MENDELIAN AND NON-MENDELIAN QUADRATIC OPERATORS 

NASIR GANIKHODJAEV, MANSOOR SABUROV, AND UYGUN JAMILOV 

Abstract. In this paper, we attempt to provide mathematical models of Mendelian and Non-Mendelian 
inheritances of the bisexual population system having Fisher's 1:1 principle. In our model, we always 
assume that distributions of the same phenotype of female and male populations are equal. We study 
the evolution of a Mendelian trait. As an application of a non-Mendelian inheritance, we construct a 
^— ^ . quadratic stochastic operator that describes transmission of ABO and Rh blood groups. 
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1. Introduction 



The study of sex allocation is often considered the most successful branch of evolutionary biology. 
Sex allocation is the allocation of resources to male and female reproductive function in sexual species 
(see ®, 0). 

In anthropology and demography, the human sex ratio is the ratio of females to males in the popu- 
lation system. Like most sexual species, the human sex ratio is approximately 1:1. Fisher's principle 
is an evolutionary model that explains why the sex ratio of most sexual species is approximately 1:1. 
It was famously outlined by R. Fisher in his book [3]. A. W. F. Edwards has remarked that Fisher's 
£ — , ■ principle is "probably the most celebrated argument in evolutionary biology" [5J. However, W.D. 
Hamilton had been introduced a model by breaking the assumptions made in Fisher's model in which 
the population system has an extraordinary sex ratio (see [5]). Sex ratios that are 1:1 are hence known 



as " Fisherian" , and those that are not 1:1 are " non-Fisherian" or "extraordinary" . 

Heredity is the passing of traits to offspring from its parent or ancestors. This is the process by which 
an offspring cell or organism acquires or becomes predisposed to the characteristics of its parent cell or 
organism. Inherited traits are controlled by genes and the complete set of genes within an organism's 
genome is called its genotype. The idea of particulate inheritance of genes can be attributed to G. 
Mendel (see [6] ) . Mendelian inheritance is a scientific description of how hereditary characteristics are 
passed from parent organisms to their offspring. In Mendelian inheritance, each parent contributes one 
of two possible alleles for a trait. In humans, eye color is an inherited characteristic and an individual 
might inherit the " brown-eye trait" from one of the parents (see [7] ) . Non-Mendelian inheritance is a 
general term that refers to any pattern of inheritance in which traits do not segregate in accordance 
with Mendel's laws. Inheritance of traits in fungi, viruses, and bacteria are all non-Mendelian. Non- 
Mendelian inheritance plays a role in several disease processes (see [8]). 

A quadratic stochastic operator (QSO) is a primary source for investigations of dynamical prop- 
erties of population genetics [9]. The fascinating applications of quadratic stochastic operators to 
population genetics were given in [TO]. It describes a distribution of a species for next generation if 
given the distribution of these species for current generation. In the paper [11], it was given a long 
self-contained exposition of the recent achievements and open problems in the theory of quadratic 
stochastic operators. 

In this paper, we are attempting to give a mathematical model of Mendelian and Non-Mendelian 
inheritances in the bisexual population systems. Our approach providing the evolution of the bisex- 
ual population system is totally different from the one given in [TO] and it gives an opportunity to 
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formulate the Mendelian inheritance in the system. We shall study the evolution of Mendelian and 
Non-Mendelian inheritances of the biosphere having Fisher's 1:1 principle. In our model, we always 
assume that " the distributions of the same phenotype of female and male populations are equal." We 
construct a quadratic stochastic operator that describe inheritance of ABO and Rh blood groups. 
This study is the continuation of the papers [12J and[13j. 

2. Preliminaries 

We recall some definitions and notions (see [12J). 

Let (A, L) be a graph without loops and multiple edges, where A is the finite set of vertexes and 
L is the set of edges of the graph and {Aj}, i = 1, ...,n— the set of maximal connected subgraphs 
(connected components) of the graph (A, L). Furthermore, let $ be a finite set, called the set of alleles 
(in problems of statistical mechanics, $— is called the range of spin). 

The function a : A — >• <£ is called a genotype (in statistical mechanics, it is called configuration). 
Denote by ft the set of all genotypes and by 5(A, $) the set of all probability distributions defined on 
the set ft. A QSO V : S(A, <J>) — > S(A, <!>) is defined as follows: for an arbitrary measure A G S(A, $), 
the measure A' = VX defined as 

(2.1) A'(<7)= Y, *WvA(a')A(<7"), Vex eft 

o-',o-"en 

where p a 'a",a = Pa"a',a > and ^ 2Vo-",o- = 1 (see [9], [10]). Here Pa'a" ,o is called a heredity 

cr£f2 

coefficient. 

In the model (|2.ip . the gender difference of species has not taken into the consideration. However, in 
many population systems, this difference plays an important role during the evolution of the population 
system. Therefore, in this paper, we are going to provide the model of the bisexual population system 
in which the gender difference has taken into the consideration. It is worth mentioning that our model 
for the bisexual population system is totally different from the model provided in |10j and it has an 
advantage to provide the mathematical model of a Mendelian inheritance. 

Let (A, V) be a graph and {Aj}, where % = 0, 1, ..., n, be sets of its maximal connected subgraphs such 
that Ao = {x }. Suppose that the population system has only female and male genders G = {/, m\. Let 
<J> be a finite set of alleles. Every species can be characterized by the following genotype a : A — » G x $ 
such that a |a : Aq — > G, and a |a\a : A \ Ao — >• $. We denote the set of all genotypes by ft. Let 
ft/ = {a G ft : <t(x°) = /} and ft m = {a € ft : a(x°) = m} be sets of all female and male genotypes, 
respectively. It is clear that ft = ft/ U ft m and ft/ n ft m = 0. For the two genotypes (parents) 
a', a" € ft, we define a set (traits of children) ft(cr',<r") such that if a' |a = °~" |a then fl(a',a") = 
and if a' |a 7^ c" |a then ft(cr', a") is some nonempty subset of ft. From the biological point of view, 
ft(er ', a") is a set of all possible genotypes (traits) of children having parents with genotypes a' and 
a". In this sense, the constrain for the set Q(a',a") means that the spices having the same gender 
would not produce anything else while the different gender spices will produce something during the 
evolution. By 5(A, G, $), we denote the set of all probability distributions of the set ft. Every element 
of the simplex 5(A, G, $) can be considered as a state of the population. 

We define the evolution of the bisexual population system as a QSO V : S(A, G, <£) — > S(A, G, <1>): 
during the evolution, every state A € £*(A, G, <3?) of the population system goes to another state 
A' = VX as follows: 

(2.2) X'(a)= Y, Pa'a",aX(a')X(a"), V a e ft 

oV"eQ 

where, p a 'a",a is heredity coefficients such that 

Pa'a''a = Pa"a',a > 
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^ / 0, if n(a',a") = 

2^ 2wv- j if n(( 7 , ,(7 ,/ )^0 

(r£fi(ff',,7") k 

^ JV<r'> = 0. 
o-£fi(<T',<r") 

Here c is a given positive constant. 

As we mentioned above that depending on possible scenarios for the evolution of the population 
system, the sex ration might or might not be 1 : 1 (see [5]). 

Therefore, in general, we suppose that the population system satisfies the following p : q law: 
during the evolution, the ratio of females and males remains the same as p : q. In other words, the 
population system has an equilibrium p : q sex ratio. 

We are going to provide some models in which one may observe a p : q law. 

We say that an evolution operator (|2.2j) has a p : q property if for any a' , a" € and (7/ S fl/, a m € 
£l m such that o-f |a\a = °m |a\a one has 

(2-3) Pa'a",a f ■ Pa'a",a m = P ■ Q, 

,„,, ^ r o if n(aV') = 

(2 - 4) I^jwvH i if iv, o 7*0, 

where, < p, q < 1 and p + q = 1. 

Biological meaning of the condition ()2.3[) in the 1 : 1 law case is that the possibilities having the 
same genotype of girls and boys are equal. 

An evolution operator (|2.2p having a p : q property is called a p : q operator. 

Let us define a hyper- simplex as follows: 

S*(A, G, $) = {Ae 5(A, G, $) : A(fy) = p, A(O m ) = g} 
One can easily check that every p : q operator has the following canonical form 
(2.5) >!{&)= 2 Yj P<T^",<rA(o J )A(o y/ ), Vaefi. 

it' en, 
CT"sn m 

Indeed, since fi = Qf U O m and ilj n fl m = 0, we have that 

A'(o-) = 5Z P CT ' CT ", CT A( ( j , )A(a ,/ ) 

o-',o-"6f2/Un m 
CT "gn m 



+ E Pa'a",aK a ') X ( a ") 

CT 'en m 
2 J] ?WvA(ct')A(<7"). 



<r"en„ 



One may check that any p : q operator satisfies a p : q /aw. 

Indeed, we want to show that the hyper-simplex <S*(A, G, <J>) is invariant under the p : q operator, 

i.e,7(S*(A,G,$))cS*(A ) G,$). 
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Let us first show that if A G 5*(A,G,$) then FA = A' G 5(A,G,$). In fact, it follows from (j2 
and (ELSD that 



J>'(<r) = 2 £ (^p^vW)A(O 

<r"e«m 

= — A(fi/) • A(fi m ) = 1. 

Now, we want to show that A' G S*(A,G, $). Indeed, for any cry G Qf,o~ m G f2 m such that 0^^ = 
0m I Ad it follows from (j2.3|) and (J2.5H that A'(<r/) : A'(cr m ) = p : q. This yields that A'(f2/) = p and 
A'(Q m ) = g. 

In the sequel, we suppose that the p : q operator acts on the hyper-simplex S*(A,G, <&), i.e., V : 
5* (A, G, $) — » S'*(A, G, $). Moreover, for the sake of simplicity, we shall consider the population 
system having a 1 : 1 law (or Fisher's principle), i.e., p = q = -^. 

3. Mendelian and Non-Mendelian QSO 

In this section, we are going to provide Mendelian and Non-Mendelian models of the bisexual 
population system having 1 : 1 law (or Fisher's principle). 

Mendel's law of the heredity can be summarized in two laws: Law of segregation and Law of 
independent assortment. Law of segregation states that when any individual produces gametes, the 
copies of a gene separate so that each gamete receives only one copy allele. Law of independent 
assortment states that alleles of different genes assort independently of one another during gamete 
formation. 

Now, we want to describe Mendelian and Non-Mendelian models of the bisexual population system. 

Let {Aj}, where i = 0,1,..., n, be a set of maximal connected subgraphs of the graph (A,L). For 
the two genotypes (parents) a', a" G ft, we define a set (traits of children) Qm(o~', <t") as follows: 
if a' |a = a" |a then S]m((t',(j") = and if a' \\ ^ a" |a then Qm(&' \ a ") = {o" G fi : u\\ i = 
o-'lAi or a\ Ai = a"\ Ai }. 

Definition 3.1. A 1 : 1 quadratic stochastic operator (12.5P is called Mendelian, if p a 'c" ,a = 0, for all 

o~ £ Qm(°~' ',?")■ A 1:1 quadratic stochastic operator (12. 5p is called Non-Mendelian if Pa'a",^ 7^ 
for some <7o G" Qm(o~' ,o-"). 

One of the natural way to construct Mendelian QSO as follows. Let S*(A, G, <£) = {A G S(A, G, <J>)A(f2/) 
A(J7 m ) = i} be a hyper-simplex. 

Let ^o £ S*(A, G, $) be a fixed measure such that 

(3.1) Vo(o-f) = Mo(o- m ), 

for any oy G fi/,<J m G U m in which ct/|a\a = O" m | A \A - 

The heredity coefficients {p a '<r",a} of the Mendelian QSO defined as follows 



2^0(0-) 

(3-2) Pa>a",a 



ifo-en M {o-',a")^ 



if a^n M (a',a")^ 

if n M (cr',a") = $ 



for any a', a" G ft. 

Then, one can easily check that 



Pa'a",a f '■ Pa'a",a m ~ 1 : 1 



^PcrV'V 



if ft(<7>") = 

2 if fi(<r',<T")^0, 
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for any a', a" G ft and aj G ft,f,o~ m G ft m such that 0/|a\a o = a m\A\A - 

This means that the 1 : 1 QSO with (13.21) is Mendelian. 

It is worth mentioning that Mendelian QSO defined by (J3.2I) is convenient to study dynamics of some 
bisexual population systems. In the next section we shall consider some applications of Mendelian 
QSO. 

4. Mendelian QSO on the simplest graph 

4.1. A Mendelian Trait. In this section we are going to discuss the mathematical model of a 
Mendelian trait. A Mendelian trait is one that is controlled by a single locus and shows a simple 
Mendelian inheritance pattern. These traits include PTC (phenylthiocarbamide) tasting, hairline 
shape, tongue rolling, earlobe attachment, hand clasping etc. (see [?]). 

Let (A,L) be a graph such that Ao = {0},Ai = {1},A = Ao U Ai and L = 0. Suppose that 
a child receives dominant and recessive alleles from parents, i.e., $2 = {^4, «}• We then have that 
O = {01, 02,0-3,0-4}, ft f = {0-1,02}, and ft m = {03,04} where 

o-i = {f,A), 2 = (/,a), a 3 = (m,A), 4 = (m,a). 

For any a' , 0" G ft one has that 



CT/ 



if0 / | Ao =o' / | Ao 

{a', a"} ifa'\ Ao ^a"\ Ao ,a'\ Al =a"\ Al 

ft if a' |a / 0" UoX |ax/ 0" |aj 



Let S*(A,G,$ 2 ) = {A G 5(A, G, 3> 2 ) : A(Q/) = A(fi m ) = ^} be a hyper-simplex. In other words, 
A G 5* (A, G, $2) means that A(oj) > 0, i = T~i such that A(oi) + A(o 2 ) = A (03) + A(o 4 ) = |. 

Now, we choose a measure //o G 5* (A, G, $2) satisfying the condition (|3.ip . i.e., /Uo(oj) = «i, i = 1, 4 
such that or = 03 = a and 02 = 04 = | — a where a G (0, |). 

We then may define heredity coefficients 0v<7" iO -} of the Mendelian 1 : 1 QSO by the formula (|3.2p . 

The evolution of the Mendelian trait of the bisexual population system can be given as follows 
A' = VX 

A'(oi) = 2A(0i)A(0 3 )+4aA(0 2 )A(0 3 ) + 

+4aA(oi)A(o4) 
A' (o 2 ) = 2(1 - 2a)A(oi)A(0 4 ) + 2(1 - 2a)A(0 2 )A(0 3 ) 

+2A(o 2 )A(o 4 ) 
(4.1) A'(o 3 ) = 2A(0i)A(0 3 )+4qA(0i)A(o 4 ) + 

+4aA(0 2 )A(0 3 ) 
A'(o 4 ) = 2(l-2a)A(0i)A(0 4 ) + 2(l-2a)A(o 2 )A(o 3 ) 

+2A(o 2 )A(o 4 ). 



Let 

(4.2) S 2n - 1 = \ x G R 2n : ^ x i = !. x i> °> Vi = V*n 



■in 



i=i 



(4.3) S*- 1 = J s G 5*- 1 : 5> = £*i = ~ 

[ i=l i=n+l J 

be a (2n — 1)— dimensional simplex and a hyper-simplex, respectively. 



(4.4) V a 
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If we denote by A(o"j) = Xi for all i = 1, 4 then the evolution operator V a : S 3 — >• 5 3 defined by (|4.ip 
has following form 

x' x = 2x1X3 + 4ax2X3 + 4axiX4, 

x' 2 = 2(1 — 2a)xiXi + 2(1 — 2a)x 2 X3 + 2x2^4, 

x' 3 = 2x1x3 + 4ax2X3 + 4axiX4, 

x 4 = 2(1 — 2a)xiX4 + 2(1 — 2a)x2X3 + 2x2X4. 

where a € (0, J?). 

It is clear that for any x € S 3 one holds (V a x)i = (V a x)3, (V a x)2 = (V a x)i. 
Let us consider the following quadratic function f a : [0, |] — > [0, ^] 

(4.5) / Q (x) = 2(1 - 4a)x 2 + 4ax 

One can easily check the following assertions: 
(i) The function (|4.5p has two fixed points Fix(f a ) = {0, ^} ; 



(ii) If a = j then f a (x 



x; 



(hi) If < a < \ then the trajectory of the function (|4.5p converges to the fixed point x* = for 
any initial point xo £ (0, 5). If 3 < a < 5 then the trajectory of the function (J4.5I) converges 
to the fixed point x* = ^ for any initial point xo € (0, ^). 

By means of the properties of the quadratic function (|4.5|) one can get the following properties of 
the operator V a : S 3 — >• S 3 given by f)4.4j) : 

(i) The operator flO|) has two fixed points Fix(V a ) = {(0, ±,0, ~) , (|,0, |,0)}; 
(ii) If a = -g then the operator (14. 4p is an identity operator; 

(iii) If < a < \ then the trajectory of the operator (14. 4p converges to the fixed point x* = 
(0, 5, 0, 5) for any initial point x° S 5 3 . If \ < a < 5 then the trajectory of the operator (|4.4|) 
converges to the fixed point x* = Q, 0, i, 0) for any initial point x° G S* 3 . 

Based on this study, we may conclude that, in the bisexual population system, the Mendelian trait 
which is dominating in numbers at the initial state would be spread towards the system in the future. 

4.2. A Mendelian Inheritance with Multiple Alleles. In this section, we are going to consider 
a Mendelian inheritance with multiple alleles. 

Let (A, L) be a graph such that Ao = {0}, Ai = {1}, A = Ao U Ai and L = 0. Suppose that a child 
receives 4 types of alleles from parents, i.e., $4 = {1, 2, 3, 4}. 
We then have that 

^ = {c/ J i,Cr/ ) 2,C"/,3)0'/,4)<7'm,l,<7'm.,2,(7'm,3,<7'm,4}, 

n f = Wf,l,Vf,2,0-f, 3 ,0-f A }, 

^m = Wm,l, 0"m,2> O m $, CTm^}, 

where 

<T ft i = (/,1), tr ft2 = (/,2), a /)3 = (/,3), a/,4 = (/,4), 
o"m,i = (m,l), <7 m ,2 = (m,2), a mj3 = (m,3), cr mj4 = (m,4). 
For any a', <r" E £1 one has that 

'0, if a' | Ao = cr" | Ao 

n^, if a' |a t^ cr" |a ,ct' |a x = c" Ui 

^cr,<5, if Cr' |a 7^ cr" Iao' "' Ui/ cr " Ui • 

where, CT = {a', a"}, £l ffi $ = {a' ,a" ,5',5"}, and 

|a = °- |Ao> d |Ai= cr |Ai 
C 1 |a = cr |a , |Ai= cr |a x • 






L M 
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Let S*(A, G, $4) = {A G S(A,G, $4) : A(S1/) = A(O m ) = i} be a hyper-simplex. In other words, 
A G 5* (A, G, $4) means that A(cj/ j j), A(cr mi j) > 0, i = 1, 4 such that 

4 4 1 



j=l 



j=l 



Now, we choose a measure fio G S^A, G, $4) satisfying the condition ()3.ip . i.e., ^o(<7/,i) = ^oi&mj) = 
a,, i = 1,4 such that a\ + . . . + a 4 = ^ and a, > for any i = 1,4. We then may define heredity 
coefficients {jVo-'v} of the Mendelian 1 : 1 QSO by the formula (|3.2p . 

The evolution of the Mendelian inheritance with 4 alleles of the bisexual population system can be 
given as follows A' = WX 

2m 



A'(o-/ ;i ) = ^ n .± n . (Hvf^HVmj) + A(cJ/ J )A(cr mi i)) 

i 



Oli + Otj 



2a 



A'(<T m ,i) = ^2 . _■_ " ( A ( <T /,j) A ( Cr "i.i) + ^(^/jO^^m.i)) 



i=i 



a* + a 



where, i = 1,4. 

If we denote by A(<7/,i) = x, and A(<7 m ,i) = Xi+4 for all i = 1,4 then the evolution operator 
VF : S* 7 — )• S 7 defined by (|4.6|) has following form 



x 'i = Yl 



la; 



an + ay 



{XiXj+4. + XjXi+4) 



i = l,4, 



'i+4 



£ 



2a,; 



«i + a,- 



(XjXj+4 + XjXj +4 ) , 2 = 1,4. 



It is clear that for any x £ S 7 one holds (Wx)i = (Wx)^^ for all i = 1,4. 
Therefore, it is enough to study the following operator 



(4.6) 



W /: 



*i = *i(2*i + ^^ + 5^x3 + ^x4), 

4 = X 2 I 2x 2 + ^xi + ^X 3 + ^X 4 ) , 
4 = X 3 (2X 3 + ^X, + ^X 2 + ^£4 ] , 



01+03' 



02+03' 



03 +04' 



4 = ^4 2X 4 + -^-X 1 + -^X 2 + 



01+0:4 



02+04' 



4o4 
03+04 



•t'3 



on the domain x\ + . . . + X4 = ^ an d x« > for all i = 1, 4. One may check that the operator (14.6 
has the following form 



(4.7) 



2(01-02) , 

01+02 



A = xi ( 1 + T^^ 2 + ^^^ xa + -^;^:' x 



2(01-03) , 
01+03 



2(01-04) , 
01+04 



2(02-01) , 

Ol+Q2 



x' = X 2 1 + ~ y ~\ " X/ Xi + ~ v ~" ~ JJ X 3 + ~ v "" " X4 



2(02—03) , 
02+03 



2(02-04) , 
02+04 



T > - Tq I 1 + 2 ( Q 3-Ol , 2(03-02) , 2(03-04) 



2(04-01) 



2(04-02) 



2(04-03), 



x 4 = x 4 i 1 + v Ql+a 7^ i + T 2 +o7 / X2 + WW 3 * 1 • 
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Let yi = 1%i and y\ = 2x[ for i = 1,4. Then the operator (|4.7p is quadratic stochastic Volterra 
operator defined on the simplex S 3 : 

(4.8) ^,( 1+ £^_^), i = - A . 

Quadratic stochastic Volterra operators defined on the finite dimensional simplex were studied 
very well in [T2]-[21J. Based on that study, we may get the same conclusion which was done for the 
Mendelian trait, i.e., in the bisexual population system, the dominant Mendelian trait of multiple 
alleles would be spread towards the system in the future. 

5. Construction of Non-Mendelian QSO on a finite graph 

In this section, we are aiming to study a special class of non-Mendelian QSO. This study will be 
applied to the transmission of Human blood group systems and Rhesus factors (see Appendix). 

Let {A,}, where i = 0, 1, ..., I, be a set of maximal connected subgraphs of the graph (A, L) such 
that Ao = {x }. Let <J> be a finite set of alleles and O be a set of all genotypes. For the two genotypes 
(parents) a', a" G CI, we define a set (traits of children) Q(o~',o~") as follows: 

(*1\ 0(„' „"\ J ' if G ' l A o =cr " lAo 

I 5 - 1 ) S2(ct,ct J = < . , . „ . 

\Ct, if a' Uo^cr |a • 

In other words, during the evolution of the bisexual population system, a child can receive any type 
of traits (not only parent's traits). 

Now, we are going to provide a construction of a non-Mendelian 1 : 1 QSO defined by a non- 
Mendelian inheritance (|5.ip . 

Let S* (A, G, $) = {AG 5(A, G, $) : X(Cl f ) = X(fl m ) = §} be a hyper-simplex. 

For given o~',cr", we choose a positive measure n a ',a" G S*(A,G, $) such that 

(5.2) /icr',a"(°"/) = M<r',er"( cr m), 

for any Of G Clf, a m G Vl m in which o/|a\a = cr m | A \ Ao . 

We define the heredity coefficients {p a 'a",a} of the non-Mendelian 1 : 1 QSO as follows 

(5 - 3) *W> = j ifT2(a>") = 0. 

Then the transmission V : 5"* (A, G, <&) — > 5* (A, G, $) of Rh blood groups can be written as follows 
A' = FA 

(5.4) A'(a)=4 ]T ^^(^A^OAK), 

<r"en m 

for all a G CI. 

Due to the condition (|5.2|) . one can see that A'(o"/) = A'(cr m ) for any cry G Qf,a m G O m in which 

°"/Ia\a = °"m|A\A - Suppose that 0/ = jo/ , ■••> °/ } and ^™ = |o"L , -,o"m J where af\ A \ Ao = 

&m |a\a for a11 i = l,n. 

Let Pi,-fc = 2/i (») o) (crl ) =2// w (j) (jL 1 , for all i,j, k = l,n. Since 

n , 

o-en k=i v 

= 1 
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n 

we have that Y^ Pij,k = 1 an d P%j,k > for alH, j, k = 1, n. 

fc=i 

We denote by A(cn ) = Xk, A(o"m ) = x n +k, VA; = 1, n. Then the non-Mendelian inheritance evolu- 
tion (15.4p can be written as follows 

n 
x k = ^ 2_/ Pij,k%i%n+j 7 

ra 
^n+fc = ^ 2^ Pij,kXiXn+j- 



for all fc = 1, n. It is clear that x' k = x' n+k . 

Let us define y^ = 2x^ and y' k = 2x' k . Then, we may study the following QSO V : S n ~ 1 — > S n ^ 1 
instead of the previous operator 



(5.5) y' k = V Pij^kViVj, k = l,n. 



n 

■ -E 



where, X] Pij.A: = 1 and Py,fc = Pji,k > ° for a11 hj, k = l,n. 
fc=i 

The regularity of the QSO (|5.5p was studied in [18]. Namely, it was proven that ifpij,k > 2k ^ or a ^ 
i, j, k = 1, n iften <5<50 ^ a « a unique fixed point and all trajectory converge to this fixed point. 

The biological interpretation of this result is that, in the bisexual population system, if initial 
distributions of all possible traits in the system are greater than the reciprocal of the total number of 
traits then eventually these distributions become stable and all types of traits will be preserved during 
the evolution. 
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Appendix A. Models of Transmission of ABO and Rh Blood Groups 

In this appendix, we construct QSO which describes inheritance of ABO and Rh blood groups. To 
the best of our knowledge, this is the first attempt to describe a mathematical model of the inheritance 
of ABO and Rh blood groups. In our model, we shall not take into consideration any internal or 
external factors of the society such as migration or immigration. We are considering an ideal society. 
We could not give any comparison between our outcomes and others because it is the first study on 
this topic. At the same time, we are not asking any credit to be true our outcomes. Of course, our 
model is far from being the best model. However, it is the first attempt in this direction and it can 
be modified and improved for any complex society. 

A.l. Transmission of Rh Blood Groups. In this section, we construct QSO which describes 
inheritance of Rh blood groups. 

It is known that there are two types of Rhesus factors: Rh positive and Rh negative, i.e., $2 = 
{+, — }. Moreover, a gender of human does not have any influence on the types of the Rhesus 
factor. Therefore, we can consider a gender of human and Rhesus factors as vertexes of discon- 
nected graph (A,L), where A = {0},Ai = {1},A = A U Ai and L = 0. We then have that 
H = {er /i+ ,o7_,o- mi+ ,o- m _}, Qf = {af t+ ,o-f-}, and Q m = {a mt+ ,a m -} where 

o-f,+ = (/,+), o-/_ = (/,-), 
a m>+ = (m, +), <7 m = (m, -). 
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Now, our main task to calculate heredity coefficients, for example p a , +1 <r m + (cr m7 _). In order to do 
so, we randomly took around 10,000 parents from Malaysia and investigated the transmission of Rh 
and ABO blood groups to their children. In our database, \pa f + ,a m , + {&m,~) means a frequency (or 
von-Mises probability) of having a son with Rh negative from both parents having Rh positive. More 

precisely, \p af + ,u m + (o~m,-) = a/7/ , + m +\ > where N(f + ,m + ) is a total number of children of parents 
both having Rh positive, and N m - (f + ,m + ) is a total number of sons with Rh negative in N(f + ,m + ). 
Similarly, we can calculate other heredity coefficients. 
For the sake of simplicity, we use the following notations 

cri = (m,+), cr 2 = (m, -), 

<?3 = (/,+), *4 = (/,-)• 

We define the all heredity coefficients p a 'c",a'- 

Pvi<Ts,*i = 0-4925 p ai a 3 ,a 2 = 0.0075 p ai a 3 ,* a = 0.4925 
P<riff 3 ,<T4 = 0.0075 Vaxo^ax = 0.3273 Pa^^ai = 0.1727 
P<riff4,<r 3 = 0.3273 p<r 1(T4j<r4 = 0.1727 2W 3 ,<ti = 0.3230 

2W 3 ,<T2 = 0.1770 P ff2(T 3 )(T 3 = 0.3230 Pa 2 a- A ,a A = 0.1770 
Po-2CT4,(Ti = 0.05 P(T2CT4,iJ2 — 0.45 Pcr2Cr4,(T3 = 0.05 

Pf72fT4,CT4 U.40 

Remark A.l. The heredity coefficients definitely depend on the region or place where we are collecting 
statistics. Consequently, the outcomes of distributions of Rh blood groups would be different from, place 
to place. In our statistics, we have randomly chosen 10,000 parents from Malaysia. Moreover, the 
size of statistics might influence to the outcomes. The outcomes become more accurate and precise as 
much as large the statistics. 

Due to (|5.4p , the transmission of Rh blood groups has the following form 

x[ = 1.9699xix 3 + 1.3094xix 4 + 1.2919x 2 x 3 + 0.2x 2 x 4 , 
, x' 2 = 0.0301xix 3 + 0.6906xix 4 + 0.7081x 2 x 3 + 1.8x 2 x 4 , 

^ '' x 3 = 1.9699xix 3 + 1.3094xix 4 + 1.2919x 2 x 3 + 0.2x 2 x 4 , 

x 4 = 0.0301xix 3 + 0.6906xix 4 + 0.7081x 2 x 3 + 1.8x 2 x 4 . 

It is clear that x\ = x' 3 and x 2 = x 4 . Let y\ = 2x\ and y 2 = 2x 2 . As we showed in the previous 
section, it is enough to study the following QSO V : S 1 — > S 1 defined on 1-dimensional simplex S 1 



(A - 2) V • I V[ = (1 - a)yl + 2(1 - b 



b)ym + (1 - c)y 2 . 

where, a = 0.9849, b = 0.6503, c = 0.1. 

The dynamics of this operator was studied in [22]. Let A = 4(1 — a)c + (1 — 2b) 2 . In our case, 
one can see that < A < 4. Then the operator (|A.2|) has a unique attracting fixed point and all its 
trajectories converge to that fixed point. Moreover, a numerical calculation shows that, a unique fixed 
point is approximately equal to (0.95, 0.05) and it is attracting. All trajectory of the operator (|A.1|) 
converge to that fixed point. 

Based on this study, we may state that the transmission of Rh blood groups in Malaysia will be 
eventually stable and Rh blood groups would be distributed as around 95% Rh positive and around 5% 
Rh negative. 

A. 2. A Transmission of ABO Blood Groups. In this section, we construct QSO which describes 
inheritance of ABO blood groups. 

It is known that there are four types of ABO blood groups: A, B, AB, O, i.e., $4 = {A, B, AB, O}. 

Moreover, a gender of human does not have any influence on the types of ABO blood groups. There- 
fore, we can consider a gender of human and ABO blood groups as vertexes of disconnected graph 
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(A, L), where A = {0}, Ai = {1}, A = A U Ai and L = 0. We then have that 

fl = {(Ti,(T2,(T3,(T4,iT5,Cr 6 ,(77,Cr 8 }, 

where 

ai = (/, A), a 2 = (/,B), a 3 = (/, AB), a 4 = (/,0), 
a 5 = (m,A), a 6 = (m,B), a 7 = (m,AB), a 8 = (m, O). 

Now, our main task to calculate heredity coefficients of ABO blood groups, for example fx fff A)(Tm A (<7 m ,A) 
In our database, fi af A)CTm A (cr m ,A) means a frequency (or von-Mises probability) of having a son 
with A blood group from parents both having A blood groups. More precisely, [i a , A ,o- m A {<J m ,A) = 

TiICF AI / All ' w here iV((/, A), (m, A)) is a total number of children of parents both having A blood 
groups, and N m) x{{f^ A), (m, A)) is a total number of sons with A blood group in N((f, A), (m, A)). 
Similarly, we can calculate other heredity coefficients. 

Here are the distributions of ABO blood groups in Malaysia: (for the sake of simplicity, we use 
fiij t k instead of ^ aiaj (cr k )) 



fjt 15jl = 0.4533 fi 16tl = 0.0865 017,1 = 0.2218 

015,2 = 0.0063 |Ui6,2 = 0.2197 0i 7 , 2 = 0.0940 

/i 15)3 = 0.0038 |Ui6,3 = 0.1661 017,3 = 0.1767 

/xi5, 4 = 0.0366 /ii 6)4 = 0.0277 017,4 = 0.0075 

fji 1Si5 = 0.4533 ^16,5 = 0.0865 017,5 = 0.2218 

/xi5 >6 = 0.0063 0i 6)6 = 0.2197 0i 7 , 6 = 0.0940 

fi 15}7 = 0.0038 /ii 6)7 = 0.1661 017,7 = 0.1767 

W58 = 0.0366 0i 6>8 = 0.0277 0i 7;8 = 0.0075 

^25,1 = 0.1750 26 ,i = 0.0060 27 ,i = 0.0192 

025.2 = 0.1833 025,2 = 0.4653 027,2 = 0.3846 

025.3 = 0.0983 ^26,3 = 0.0060 //27,3 = 0.0769 
1X25,4 = 0.0433 ^26,4 = 0.0227 027,4 = 0.0192 
025,5 = 0.1750 26 ,5 = 0.0060 27 ,5 = 0.0192 
^25,6 = 0.1833 ^6,6 = 0.4653 ^7,6 = 0.3846 
^ 257 = 0.0983 ^26,7 = 0.0060 27 , 7 = 0.0769 
^25,8 = 0.0433 ^6,8 = 0.0227 ^7,8 = 0.0192 
M351 = 0.2525 /^i = 0.0714 37 ,i = 0.0522 

035.2 = 0.0707 03 6i 2 = O.2662 37 , 2 = 0.0373 

035.3 = 0.1667 036,3 = 0.1299 037,3 = 0.4030 

035 4 = 0.0101 036 4 = 0.0325 037 4 = 0.0075 
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//35,5 = 0.2525 ^36,5 = 0.0714 ^7,5 = 0.0522 

^35,6 = 0.0707 ^6,6 = 0.2662 ^7,6 = 0.0373 

M357 = 0.1667 ^36,7 = 0.1299 ^37,7 = 0.4030 

M35) 8 = 0.0101 ^36,8 = 0.0325 ^37,8 = 0.0075 

/i 4 5,i = 0.2730 ^46,1 = 0.0085 ^47,1 = 0.0709 

^45,2 = 0.0117 ^46,2 = 0.2888 ^47,2 = 0.0613 

M453 = 0.0059 ^46,3 = 0.0049 ^47,3 = 0.1142 

/M5,4 = 0.2094 fx 46A = 0.1977 ^47,4 = 0.2536 

/i 4 5,5 = 0.2730 ^46,5 = 0.0085 ^47,5 = 0.0709 

^45,6 = 0.0117 ^46,6 = 0.2888 ^47,6 = 0.0613 

^45,7 = 0.0059 /u 46 ,7 = 0.0049 ^47,7 = 0.1142 

/i 4 5,8 = 0.2094 ^ 46j8 = 0.1977 ^47,8 = 0.2536 

/ii 8 ,i = 0.2146 /u 28 ,i = 0.0066 /z 38 ,i = 0.1070 

Mi8,2 = 0.0022 ^28,2 = 0.2325 ^38,2 = 0.1209 

/ii8, 3 = 0.0045 ^28,3 = 0.0066 ^38,3 = 0.0651 

/ii 8 ,4 = 0.2787 ^28,4 = 0.2542 ^ 38 , 4 = 0.2070 

/il8,5 = 0.2146 ^28,5 = 0.0066 ^38,5 = 0.1070 

^8, 6 = 0.0022 ^8,6 = 0.2325 ^8,6 = 0.1209 

/ii8, 7 = 0.0045 /u 2 8,7 = 0.0066 ^ 38 , 7 = 0.0651 

^8,8 = 0.2787 ^8,8 = 0.2542 ^8,8 = 0.2070 

M481 = 0.0058 /u 48 ,2 = 0.0035 ^ 4 8,3 = 0.0034 

M484 = 0.4873 /u 48 ,5 = 0.0058 ^ 4 8,6 = 0.0035 

HMj = 0.0034 ^ 48i8 = 0.4873 

Remark A. 2. The heredity coefficients of ABO blood groups definitely depend on the region or place 
where we are collecting statistics. Consequently, the outcomes of distributions ofof ABO blood groups 
would be different from place to place. In our statistics, we have randomly chosen 10,000 parents from 
Malaysia. Moreover, the size of statistics might influence to the outcomes. The outcomes become more 
accurate and precise as much as large the statistics. 

Due to (j5.4|) . the transmission V : S 7 — >• S 7 of ABO blood groups has the following form 

x\ = 1.8131xix 5 + 0.3460xix 6 + 0.8872xiX7 + 0.858xia; 8 

+ 0.7x2x5 + 0.0239x 2 x 6 + 0.0769x 2 x 7 + 0.0265x 2 x 8 

+ I.OIOIX3X5 + 0.2857x 3 x 6 + 0.209x 3 x 7 + 0.427x 3 x 8 

+ 1.092x 4 x 5 + 0.0339x 4 x 6 + 0.2837x 4 x 7 + 0.023x 4 x 8 

x' 2 = 0.0253xix 5 + 0.8789xix 6 + 0.3759xix 7 + 0.009xix 8 

+ 0.7333x 2 x 5 + 1.8612x 2 x 6 + 1.5385x 2 x 7 + 0.930x 2 x 8 

+ 0.2828x 3 x 5 + 1.0649x 3 x 6 + 0.1493x 3 x 7 + 0.483x 3 x 8 

+ 0.047x 4 x 5 + 1.1554x 4 x 6 + 0.2452x 4 x 7 + 0.014x 4 x 8 
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0.0152xia; 5 + 0.6644xix 6 + 0.7068xix 7 + 0.018xix 8 
0.3933x 2 x 5 + 0.0239x 2 x 6 + 0.3077x 2 x 7 + 0.026x 2 x 8 
0.6667x 3 x 5 + 0.5195x 3 x 6 + 1.6119x 3 x 7 + 0.260x 3 x 8 
0.0235x 4 x 5 + 0.0198x 4 x 6 + 0.4567x 4 x 7 + 0.013x 4 x 8 

0.147xix 5 + O.lllxixe + 0.0301xix 7 + 1.115xix 8 
0.1733x 2 x 5 + 0.0909x 2 x 6 + 0.0769x 2 x 7 + 1.017x 2 x 8 
0.0404x 3 x 5 + 0.1299x 3 x 6 + 0.0299x 3 x 7 + 0.828x 3 x 8 
0.8376x 4 x 5 + 0.791x 4 x 6 + 1.0144x 4 x 7 + 1.949x 4 x 8 

1.8131xix 5 + 0.3460xix 6 + 0.8872xix 7 + 0.858xix 8 
0.7x 2 x 5 + 0.0239x 2 x 6 + 0.0769x 2 x 7 + 0.026x 2 x 8 
1.0101x 3 x 5 + 0.2857x 3 x 6 + 0.209x 3 x 7 + 0.427x 3 x 8 
1.092x 4 x 5 + 0.0339x 4 x 6 + 0.2837x 4 x 7 + 0.023x 4 x 8 

0.0253xix 5 + 0.8789xix 6 + 0.3759xix 7 + 0.009xix 8 
0.7333x 2 x 5 + 1.8612x 2 x 6 + 1.5385x 2 x 7 + 0.930x 2 x 8 
0.2828x 3 x 5 + 1.0649s 3 S6 + 0.1493x 3 x 7 + 0.483x 3 x 8 
0.047x 4 x 5 + 1.1554x 4 x 6 + 0.2452x 4 x 7 + 0.014x 4 x 8 

0.0152xix 5 + 0.6644xix 6 + 0.7068xix 7 + 0.018xix 8 
0.3933x 2 x 5 + 0.0239x 2 x 6 + 0.3077x 2 x 7 + 0.027x 2 x 8 
0.6667x 3 x 5 + 0.5195x 3 x 6 + 1.6119x 3 x 7 + 0.261x 3 x 8 
0.0235x 4 x 5 + 0.0198x 4 x 6 + 0.4567x 4 x 7 + 0.014x 4 x 8 

0.1465xix 5 + 0.1107xix 6 + 0.0301xix 7 + 1.115xiX 8 

0.1733x 2 x 5 + 0.0909x 2 x 6 + 0.0769x 2 x 7 + 1.017x 2 x 8 

0.0404x3X5 + 0.1299x 3 x 6 + 0.0299x 3 x 7 + 0.828x 3 x 8 

0.8376x 4 x 5 + 0.791x 4 x 6 + 1.0144x 4 x 7 + 1.949x 4 x 8 

This operator has a unique fixed point 

(0.042, 0.258, 0.029, 0.171, 0.042, 0.258, 0.029, 0.171) 

which is attracting and all its trajectory converge to this fixed point. 

Based on this study, we can state that the transmission of ABO blood groups in Malaysia will be 
eventually stable and ABO blood groups would be distributed as follows: around 8% from A, around 
52% from B, around 6% from AB, and around 34% from O. 
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